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Abstract. On Monday January 19, 2015 a story broke that the National 
Football League (NFL) had started an investigation into whether the New 
England Patriots deliberately deflated the footballs they used during their 
championship win over the Indianapolis Colts. Like an infectious disease, 
discussion regarding Deflategate grew rapidly on social media sites in the hours 
and days after the release of the story. However, after the Super Bowl was 
over, the scandal slowly began to dissipate and lost much of the attention it 
had originally had, as interest in the NFL wained at the completion of its 
season. We construct a simple epidemic model for the infectiousness of the 
Deflategate news story. We then use data from the social media site Twitter 
to estimate the parameters of this model using standard techniques from the 
study of inverse problems. We find that the infectiousness (as measured by 
the basic reproduction number T^o) of Deflategate rivals that of any infectious 
disease that we are aware of, and is actually more infectious than recent news 
stories of greater importance - both in terms of T^o and in terms of the average 
amount of time the average tweeter continued to tweet about the news story. 
Keywords, disease model, inverse problem, national football league, basic 
reproduction number 


1. Introduction 

The National Football League (NFL) has made American football arguably the 
most popular sport throughout the United States. The NFL was formed in 1922 
from the American Professional Football Association, and originally consisted of 18 
teams. By 1925 the league began drawing tens of thousands of fans into stadiums 
to watch games live, and by 1934 the first NFL game between the Chicago Bears 
and the Detroit Lions on Thanksgiving Day was broadcasted live on national radio - 
allowing the league’s fan base to spread nationally. The NFL’s popularity continued 
to grow when the Brooklyn Dodgers played the Philadelphia Eagles was the first 
NFL game to be televised on the National Broadcasting Company (NBC). Numerous 
other professional football leagues were formed to compete with the NFL during 
this time - although none could compete against the NFL until the All-America 
Football Conference (AFC), which later merged with the NFL in 1949 to further 
increase the size and fan base of the NFL. As the radio and television became more 
widespread across the United States, the NFL was able to increase their popularity 
as almost all NFL games were broadcasted on radio or television by 1964, allowing 
professional football to overtake professional baseball in popularity around 1965. 
As the years progressed, the popularity of the NFL continued to grow as television 
ratings surged into the millions viewers and stadiums grew to house over 100,000 
fans. Today, the fan base of the NFL well over 100 million people, with increasing 
popularity overseas [12]. 
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Despite the popularity of the NFL, its history has been plagued by numerous 
scandals; most recently, a scandal involving the footballs used by the New England 
Patriots in the 2015 AFC Championship game against the Indianapolis Colts, which 
has become known as “Deflategate”. On Monday January 19, 2015 a story broke 
that the NFL had started an investigation into whether the New England Patriots 
deliberately deflated the footballs they used during their AFC Championship win 
over the Indianapolis Colts. The NFL rules state that footballs must weigh between 
14 and 15 ounces and be inflated between 12.5 and 13.5 pounds per square inch. 
Deflated footballs would have given the Patriots an advantage over the Colts (as 
well as other opponents) by providing a better grip on the ball [8]. 

Discussion was sparked instantly on social media regarding Deflategate, and 
grew rapidly in the hours and days after the release of the story. This scandal 
gained national attention quickly as the New England Patriots had just earned a 
trip to Super Bowl XLIX against the Seattle Seahawks (another team displaying 
arguable ethics of its own in previous seasons m)- Since the Super Bowl is the 
most anticipated game in the NFL season, attention to the story was heightened. 
After the Super Bowl was over, the scandal slowly began to dissipate and lost 
much of the attention it had originally had, as interest in the NFL decreased at 
the completion of its season.This interest continued to dissipate until early May, 
2015, when the aftermath of a 243-page report by independent attorney Theodore 
V. Wells, Jr. resulted in the NFL suspending Patriots star quarterback Tom Brady 
for the first four games of the 2015 season and stripping the team of two high draft 
picks. 

The sharp rise in interest in the Deflategate scandal and then an initial slow 
dissipation of interest is similar to that of an outbreak of an infectious disease. 
When an infectious disease is first introduced to a new population, it can spread 
rapidly as a large number of the population becomes infected with the disease. The 
spread of the disease slows to more manageable rate as more of the population 
becomes immune to the disease through vaccinations or survival of victims of the 
disease, leaving less of the population susceptible to the disease. The infectiousness 
of a disease can often be captured by a dimensionless parameter called the basic 
reproduction number TZq (“R-naught”). The value of TZq informs scientics of how 
fast the disease is initially spreading by indicating the number of individuals an 
initial infected individual is likely to infect. For example, the TZq value for the 
Ebola virus was 2, meaning that if one person was infected with Ebola during the 
early stages of the epidemic, they would likely spread the virus to 2 other people 
before they recovered or passed away. Some other TZq values include 4 for the HIV 
and SARS virus, 10 for mumps and 18 for measles An TZq value of 18 is 

extremely high and means that measles was extremely infectious; which is why it 
was so catastrophic. 

Since the Deflategate scandal spread in a way similar to an infectious disease, 
we can derive and use something like an TZq value for this news event in order 
to understand how infectious this story was. To do this we must determine a 
useful medium through which information was spread, which we determined to be 
the social networking site Twitter (http://twitter.com). Twitter allows users post 
tweets, a message under 160 characters, to their Twitter page so that their followers, 
other users who will see the tweet, can read their message. The followers can then 
reply to the message and start a conversation with the owner of the tweet. They 
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can also simply retweet that tweet to their followers. A retweet is when follower 
takes a tweet they saw and posts it onto their Twitter page for their followers to 
see. The number of tweets and retweets about Deflategate can be used to determine 
the TZq value and the infectiousness of this scandal as a news story. 

In this paper we construct a simple SIR epidemic model for the infectiousness of 
the NFL’s Deflategate news story. We use standard techniques from epidemiology 
to derive the basic reproduction number TZq for this model, as well as the average 
amount of time the average individual tweets about the story. We then use Twitter 
data to estimate the parameters of this model using standard techniques from the 
study of inverse problems. We find that the infectiousness (as measured by TZq) of 
Deflategate rivals that of any infectious disease that we are aware of, and is actually 
more infectious than the story of Hillary Clinton’s announcement of her presidential 
campaign in April, 2015 - both in terms of TZq and in terms of the average amount 
of time the average tweeter continued to tweet about the news story. We also show 
that the average individual tweeted about Deflategate roughly ten times more than 
they did about the story of Freddie Gray’s death that elicited the Baltimore riots 
in April of 2015. 


2. The Model 

While the punishment of the Patriots in May, 2015 created its own news story 
regarding Deflategate, we will consider only the dynamics of the initial story, occur- 
ing in January, 2015. We assume that information concerning Deflategate is spread 
person to person like an acute infection. Since online discussion of specific topics 
is often measurable, we will focus on the subpopulation of media consumers that 
use the social media site Twitter. Time t will be measured in days since the first 
report of Deflategate on January 19th, 2015. Using formalism first studied in depth 
by Kermack and McKendrick [TH], and popularized by Brauer and Castillo-Chavez 
[5], to name just two, we categorize the individuals in the Twitter population as 
susceptible, infectious, and removed. The susceptible population S{t) (measured in 
thousands) are those individuals that regularly read and comment on sports news 
but have yet to comment on Deflategate as of time t. The infectious population I(t) 
(measured in thousands) are those that are tweeting (or retweeting) posts about 
Deflategate using the keywords #deflategate, deflate gate, deflate-gate, spygate, or 
“deflated balls.” As a simplifying assumption we consider each tweet as represent¬ 
ing a unique individual (tweeter), and thus tweets and individuals in the infectious 
class are discussed interchangably. The removed population i?(t) (measured in thou¬ 
sands) includes either those individuals who do not read or comment on sports (and 
are hence removed from this news event even from the very beginning), or individ¬ 
uals that were once tweeting about the story and have permanently stopped using 
the Deflategate keywords. Due to the short lifespan of the story (in the order of 
weeks), we assume that the immigration of new users, emigration of current uses, 
birth, and death can be neglected, yielding a constant total population size. 

The progression from the susceptible population S to the infectious population I 
to the recovered (or - maybe more appropriately named - bored) population R can 
be visualized in the traditional SIR conceptual diagram (Figure [1]). The progres¬ 
sion of the information through the population depends on many factors. One of 
the most prominent is the total number of “interactions” between susceptible and 
infectious individuals. Information about Deflategate spreads when a susceptible 
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Figure 1 . Conceptual model of information transmission and re¬ 
covery into the removed population. Black arrows show the move¬ 
ment between the S and I classes and the I and R classes. The 
fact that the level of the size of the infection population influences 
the rate at which a susceptible individual moves in to the infectious 
class is show by the dashed arrow. 


individual comes in contact with the information spread by an infectious individual 
(by reading his or her tweets) and subsequently becomes infectious (starts tweeting 
themselves). Mathematically, a reasonable measure for the number of encounters 
between susceptible and an infectious individuals, assuming homogeneous mixing, 
is given by the product SI. This is referred to as the law of mass action in the 
applied mathematics literature m- However, not every interaction of a suscep¬ 
tible person reading a tweet about Deflategate will cause a retweet or a series of 
original tweets by the newly-infectious. We use a parameter /3, the transmission 
coefficient, as a daily measure of the probability that a single susceptible reader will 
retweet using the keywords #deflategate, “deflate gate”, “deflate-gate”, “spygate” 
(a previous Patriots scandal), or “deflated balls” upon encountering such tweets in 
his or her feed. Other functional forms exist to model the rate of new infections 
(for example, the Holling Type II functional form |15jl. which we will also explore. 
Once infectious, we assume that individuals tweets about this story for an average 
of 1/7 days. The following assumptions can be collected to create the following 
SIR model 

S'it) = 

i'{t) = ps{t)m-^m ( 1 ) 

5(0) =5o, /(0)=/o, 

where 5o and Iq are parameters giving the initial population of susceptible and 
infectious individuals (measured in thousands). Since the Rif) population doesn’t 
interact with the rest of the population at any given time, we can omit it from our 
analysis. We assume that the parameters /3, 7 , 5o and Iq are all positive, which 
implies that the solution vector [Sf) If)]^ exists and is positive for alH > 0 pi] . 

The derivation of the basic reproduction number TZq for this SIR model follows 
from linearizing about the “disease-free” state [5o 0 ]"^, where the susceptible pop¬ 
ulation has yet to be exposed to the news story. TZq will represent the number of 
susceptibles that are initially “infected” by an individual already tweeting about 
the story and whether the story will spread rapidly and receive attention or will 
receive little attention causing the story to fall from the headlines and die off. 

According to the SIR modellU Sf) will always be a decreasing function. However 
for a small initial tweeting population Iq, the rate of change of S maybe small and 
therefore we assume Sf) « Sq- Consequently, the equation for If) reduces to 

I'f) « pSoIf) - ilf) = (/35o - i)If) = klf) 


(2) 
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where fc is a constant that is equal to f3So — 7- The solution to the linear equation 
([2]) is I{t) loe^*. The value of k determines whether the function increases or 
decreases for small t. The infected population will initially increase when k = 
/3So — j > 0. Or equivalently, 

,3So - 7 > 0 
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Similarly the infected population will decrease if k = PSo/j < 1. This threshold is 
the basic reproduction number TZq = PSo/"f. One can interpret this value as the 
total rate of initial infections (PSo) times the average amount of time spent infected 
( 7 “^). Thus, for fixed Sq, diseases with large P and/or small 7 will be the the most 
likely diseases to result in epidemics. 

3. Methods 

The size of the infected population during the days after Deflategate was first 
reports was gleaned from an articlej^ published in Boston Magazine. Data pre¬ 
sented in the article was compiled by the website Topsy (http://www.topsy.com), 
an analytics company and certified Twitter partner, that collects Twitter data over 
the span of 30 days. The number of tweets per day using the keywords #deflate- 
gate, “deflate gate”, “deflate-gate”, “spygate”, or ’’deflated balls” is given by the 
unfilled circles in Figure 2. 

We used this Twitter data to reverse-engineer the values of the parameters 
P, 7 , Sq and Iq in the model o by using standard methods from the study of 
inverse problems |3] , which we summarize below. 

To employ the techniques from inverse problems we require a statistical model 
to go along with the mathematical model dH). To create such a statistical model, 
we abstract our mathematical model o as in [ 2 ] to give us 

x'{t,0) = g{x(t,0),0), t€[to,tf], (3) 

x{to,0) = Xq 

where x{t,0) = I{t)]'^ and xq = [<S'o, Iq]^ are the vectors of state variables 

and initial conditions of our system, given the parameter vector 0 = [/3, a, Sq, /q]^. 
We define as our observation process the following 

f{t,0) = Cx{t,0)=I{t), 

as we are only aware of the infectious population (i.e. the tweeters) at any given 
time t. In this case C = [0 1] is a functional over K.^ - the case where we were able 
to see the entire system [S{t) I{t)]'^ at any given time the operator C would be 
the 2 x 2 identity matrix. 

To find an estimate for the parameter vector 0 we formulate the statistical model 
Y{t) = f{t,0Q) + e{t) 

= lit,0o) + e{t), t€[to,tf], 

where Oq = [Pq, ap, 5 * 00 , /po] is considered the vector of hypothesized “true values” 
of the unknown parameters and e(t) is a random variable that represents observation 
error for the observed state variable at each time t. We assume that the error 
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function e is such that E{e(t)) = 0, Var{e{t)) = cr^ and Cov{e(t)e{s)) = a'^S(t — s) 
for all t,s £ [to,tf]. 

Since our data is collected at discrete time points, we have to write our above 
statistical model in discrete terms. Given data Ii, In taken at time points 
to < ti < t 2 < < ... < tn < tf we write our observation process as 

j = l,2,...,n, 

with the discrete statistical model as 

X?' “ /(^ji^o) + ^{tj) = I{tj,6) + e{tj) 
and a realization of the discrete statistical model as 


h - + G - + G- 

Knowing the value 9 and solving the system of differential equations (ED with these 
parameter values is known as the forward problem. Alternatively, having a set of 
data /i, I 2 , In and estimating 9o is known as solving an inverse problem. In this 
paper we perform the latter. 

While several methods exist to solve inverse problems [5], we will use ordinary 
least squares. In addition to the assumptions above, we assume that realizations 
of e{t) at particular time points are independent and identically distributed normal 
random variables. In this case, one can show that the parameter vector 60 that 
maximizes the likelihood function 

C{9o) = P{Ij = I{tj,e),j = 1 , 2 ,..., n\e = Oo) 
is given by the parameter vector that minimizes the least-squares functional 

n 

CC0O) = Y.{i, -i{t„e))\ 

i=i 


In other words 

Oo = argmin^gQ£/:(0), (4) 

where 0 is the set of admissible values for the parameter vector 9. For the model 
in this paper 


0 = (0,oo) X (0,oo) X (0,3.02 x 10®) x (0,3.02 x 10®), 


where 3.02 x 10® are the estimated number of monthly Twitter users as of 5/14/2015 

PD- 

Once an estimate for parameter vector 6*0 is found, the next question is in regards 
to the uncertainty in the estimate. To generate standard errors for each of the 
parameters we create 1x4 sensitivity matrix 


Dj{9) 


dl{tj,9) dl{tj,9) dl{tj,9) dl{tj,9) 
d9i 802 302 802 


from which we create the 4x4 convariance matrix 


E"(0o) = (a2)-i I Y,Df;{0o)D,{0o) 


It follows that the standard error in the kth. component of the parameter vector 0o 
can be approximated by the square root of the fc, fcth element of the matrix I]”( 0 o) 
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Parameter 

lower bound 

estimate 

upper bound 

/3 

0.0130 

0.0153 

0.0177 

7 

0.2830 

0.3643 

0.4456 

^0 

139.0034 

156.6120 

174.2207 

.^0 

0.6626 

2.2726 

3.8827 


Table 1. Parameter values and their respective 95% confidence intervals. 


[2]. We report parameter values in an interval with the lower and upper bounds 
being two standard errors below and above the estimated value, respectively. 

All programming was performed in R |20j . with code available upon request. 
To solve this inverse problem we used a fourth-order Runge-Kutta [53] to solve 
the forward problem ([3]) for each set of parameters 0 in 0, and used the built-in 
optimization algorithm “optim” in R [20] to solve the minimization problem ([¥|) for 
Oo- 


4. Results 

For the data in Figure 2 a parameter vector that minimized the value of CC is 
given by 

00 = [Po 70 5*00 loof = [0.0153 0.3643 156.6120 2.2723]^. 

Both and S {t) and I(t) components of the solution to (|T|) subject to these parameter 
values is displayed in Table 1, along with the Twitter data used to solve the inverse 
problem. 95-percent confidence intervals for the parameters are given in Table 1. 
Since all of these confidence intervals exclude zero, we can be fairly certain that 
each of the parameters in the model are significantly different than zero, and thus 
necessary to include. While there were other parameter values that were (local) 
minimums of CC in the space 0, the fit in Figure 2 suggests that 6 o is a reasonable 
estimate for the parameter vector 9. 

The parameter values in 9q are able to give us quantitative information regarding 
the infectiousness of the Defiategate news story. For example, the estimate for /3 
suggests that between 1.3 and 1.7 percent of susceptible tweeters’ views of Defiate¬ 
gate tweets will result in the immediate conversion of said tweeter into someone 
tweeting about the Defiategate story. The estimate for 7 suggests that the av¬ 
erage Defiategate tweeter spends between 2.24 and 3.53 days tweeting about the 
Defiategate story before becoming bored with the story. 

Since our units for S and I are in terms of thousands of tweeters, our value 
for the basic reproduction number TZq = actually corresponds to the 

average number, in thousands, of new Defiategate tweeters elicited by the aver¬ 
age initial thousand Defiategate tweeters. Our estimates in do therefore suggest a 
basic reproduction number between 4.05 and 10.91 thousands of new Defiategate 
tweeters per initial thousand Defiategate tweeters, a number rivalling some of the 
aforementioned epidemics in history. 

When viewing S and I together on the same set of axes we see that by the sixth 
day of the news story the entire population has either become a Defiategate tweeter 
themselves or was never going to tweet about Defiategate in the first place. On the 
other hand, the values of / continue to be well above zero almost two weeks into 










ERIC EAGER*’'^, MEGAN EBERLE"^, AND JAMES PEIRGE*’'^ 



t 

Figure 2. The number (in thousands) of Twitter users suscepti¬ 
ble, S{t), and tweeting about Deflategate, I{t), subject to param¬ 
eter vector 9 solving the inverse problem using Topsy.com data 
(unfilled circles), t = 0 corresponds to January 17, 2015. 


the news story. These two pieces of information suggest that this news story - and 
possibly NFL news stories in general - are quite infectious and have a relatively 
high amount of staying power, which coincides with our initial intuition |12] . 

We also explored the possibility of another functional form for the probability of 
an “infection”, using the Holling type II [15] functional form 




psi 

l + hS 


in place of traditional mass action PSI, modeling the chance that too many sus¬ 
ceptible individuals can cause a saturation effect in dynamics of the news story. 
Here, h is known as the “handling time” in population ecology, or the amount of 
time it takes for the average “infection” to occur. In the context of a news story, 
this could be thought of as the amount of time it takes for an individual to realize 
the story was important enough to tweet about. When solving the inverse prob¬ 
lem in this case we found that the 95 percent confidence interval for h was equal 
to (—0.0005,0.0034), which includes h = 0, suggesting that the infection medium 
(social media) is such that information transfer moves sufficiently fast so as not 
to require any handling time, or that this NFL story was considered “important” 
roughly immediately upon the converted tweeter’s read of the infectious tweet, 
eliciting a subsequent tweet shortly thereafter. 
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5. Discussion 

In this paper we created, analyzed and fit a model for the infectiousness of 
the NFL’s initial Deflategate story using a standard SIR epidemiological model 
and data from the population social networking site Twitter. We found that this 
standard SIR model fit our data quite well (see Figure [2]), suggesting that the as¬ 
sumptions inherent in using such a model are reasonable. In fact, there have been 
some studies (for example, [19] and i) suggesting that mass-action assumptions 
may be improperly applied in the study of actual epidemics when the underlying 
populations are a) too small b) not homogeneous in space enough to warrant such 
a simple transmission probability or c) too crowded so as to saturate infectious¬ 
ness when pushed beyond a certain population size. Populations of Twitter users, 
however, consist of many individuals on one webpage unimpeded physical limited, 
alleviating the aforementioned issues. Thus, the simplest SIR model may, in many 
ways, be a better initial model for some instances of information moving through a 
social networking site like Twitter than it is for an epidemic moving through a real 
population. 

We used the parameters from the study of this inverse problem to determine how 
popular and persistent this news story was in terms of the composite parameter 
Rq and the parameter 7 , respectively. We found that the average group of 1000 
individuals tweeting about Deflategate were able to elicit 4.05 to 10.91 thousand 
new tweeters to tweet about Deflategate during the early stages of the news story, 
and that the average individual tweeting about the story tweeted between 2.24 and 
3.53 about the story. To put these numbers in perspective, we ran a similar anal¬ 
ysis of two different, prominant news stories that happened since the Deflategate 
story: the announcement of Hillary Clinton’s presidential campaign and the stories 
surrounding the Baltimore riots resulting from the death of Freddie Gray. Inter¬ 
estingly, the SIR model used in this paper was not able to capture the behavior of 
either story, due to the almost immediate emergence of individuals tweeting about 
the story (see Figure[3]and FigureHj). Instead of using the traditional SIR model in 
these cases we used an SIR model allowing for recruitment of new susceptibles for 
the Clinton story and used a simple exponential model in the Freddie Gray case. 
In the former we are still able to recover TZq and 1 /y and the latter only an upper 
bound for I/ 7 . We found that Hillary Clinton’s annoucement of her presidential 
bid, while having more total tweets at its peak than the Deflategate story (see Fig¬ 
ure 3), and a similar initial infectiousness {TZo G (4.95,7.69)), had far less staying 
power. The average tweeter was only tweeting about Mrs. Clinton’s annoucement 
for between 0.539 and 0.7133 days. This may be due to the fact that the election 
was still more than a year away, or that there were other candidates announcing 
their bids during the same period of time. Be that as it may, the model for the 
Baltimore riot story revealed something far more disturbing. This story, while elic¬ 
iting far more tweets than the Deflategate story, saw the average tweeter tweeting 
about the riots for a time whose upper bound is between only 0.3227 and 0.3232 
days (see Figure 4). One explanation for this short staying time would be that, 
after many instances of questionable police actions over the course of the past year, 
many people are starting to grow fatigued by such stories. However, one would 
hope for the same such fatigue over relatively meaningless NFL stories, which we 
have shown is not the case. 
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Figure 3. The number (in thousands) of Twitter users tweeting 
about Hillary Clinton’s announcement for candidacy for the 2016 
presidential race subject to the parameter vector 6 solving the in¬ 
verse problem using Topsy.com data (unfilled circles). We used an 
SIR model with a sorce term for the susceptible population (i.e. 
S'{t) = r — pS{t)I{t)) to better fit the data, t = 0 corresponds to 
April 7, 2015. 


The results in this manuscript suggest that the NFL’s popularity rivals (even 
surpasses) that of what should be some of our nation’s (and world’s) biggest news 
stories. Observations made in the wake of various news stories involving the in¬ 
tersection of a game (football) with some of our nation’s most contentious issues 
(gambling, equality, domestic abuse, child abuse and financial literacy, to name a 
few) appear to coincide with these results ([H], [1], [10], ^2\, [E], [S])- Even Noam 
Chomsky, in an exerpt from The Chomsky Reader, has pondered “why Americans 
know so much about sports but so little about world affairs”, arguing that sports 
provide the easiest way for people to concentrate each other (into fan bases) for a 
common (if inconsequential) goal or purpose |7]. This large-scale societal escape, 
he argues, is to combat the feeling shared by many that they simply don’t have 
much power to solve society’s real problems. However, the rising popularity of 
the NFL, along with the societal stories and scandals that its members (players, 
coaches, execs) create, if properly leveraged, can provide society with the means to 
see societal change through the inevitable clash between humanity and sport. 
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Figure 4. The number (in thousands) of Twitter users tweeting 
about Freddie Gray and the Baltimore Riots subject to parameter 
vector 9 solving the inverse problem using Topsy.com data (unfilled 
circles). Only the second half of the data (starting with t = 6) was 
fit to a simple exponential model due to the inability to fit an SIR 
model to all of the data, t = 0 corresponds to April 21, 2015. 
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